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Semilocal and dispersion-corrected density-functional calculations have been performed to study 
the crystal structure, equation of state, and electronic structure of metal tellurohalides with chem¬ 
ical formula MeTel where Me=Bi, An, or Pd. A comparative investigation of the results of these 
calculations is conducted, which reveals the role of van der Waals attraction. It is shown that the 
prediction of crystal structure of metal tellurohalides is systematically improved thanks to the in¬ 
clusion of van der Waals dispersion. It is found for BiTel and AuTel that the energy versus volume 
curve is anomalously flat in the vicinity of equilibrium volume and the calculated equation of state 
has an excessively steep slope in the low-pressure region, which are also fixed in the dispersion- 
corrected calculations. Analysis based on the computation of the volume and axial compressibilities 
shows that predicting the anisotropy of BiTel via the semilocal calculations yields an unrealistic 
result whereas the results of dispersion-corrected calculations agree with the experimental com¬ 
pressibility data. Our calculations render that BiTel (AuTel) is a narrow band gap semiconductor 
with Rashba-type spin-splitting at the band edges (with an indirect band gap) while PdTel is a 
metal with relatively low density of states at the Fermi level. The band gaps of BiTel and AuTel 
obtained via semilocal (dispersion-corrected) calculations are found to be greater (smaller) than 
the respective experimental values, which is against (in line with) the expected trend. Similarly, 
the Rashba parameters of BiTel are bracketed by the respective values obtained via semilocal and 
dispersion-corrected calculations, e.g., a larger value for the Rashba parameter or is obtained in 
association with the reduction of the band gap caused by modification of the crystal structure ow¬ 
ing to van der Waals attraction. Excellent agreement with the experimental Rashba parameters is 
obtained via interpolation of the calculated (semilocal and dispersion-corrected) values. 

PACS numbers: 71.20.Ps,71.70.Ej,64.30.Jk,31.70.-f 


I. INTRODUCTION 

Bismuth tellurohalide (BiTel) has recently been at¬ 
tracting a great deal of interest as a potential spin- 
tronic material due to emergence of giant Rashba-type 
spin-splittings in its band structure^ and of a pressure- 
induced topological insulating phase,albeit it has also 
been pointed out^ that the latter would probably be hin¬ 
dered by a structural phase transition. Prompted by 
the discovery^ of Rashba splittings in the band struc¬ 
ture of BiTel, density-functional calculations have in¬ 
creasingly been employed®^^^ to explore the electronic 
and crystal structure of this semiconductor. Curiously, 
although BiTel is often described to have a layered struc¬ 
ture where the triple Bi-Te-I layers are stacked along an 
axis perpendicular to the layers, the van der Waals at¬ 
tractions were not explicitly included in these studies. 
The lattice parameters of BiTel were therefore needed to 
be fixed to their experimental values^^ in most density- 
functional calculations. In order to provide a more com¬ 
plete description, we adopt here the dispersion-corrected 
density-functional (DFT-D2) approach,^® where the van 
der Waals forces are incorporated by means of a semiem- 
pirical force field. We employ the DFT-D2 calculations 
for not only BiTel but also AuTel and PdTel as far as we 
are aware of few experimental studies^^^^® on the struc¬ 
ture and properties of the latter, which have not been 
explored via first-principles methods before. 


In order to reveal the effect of the dispersion correc¬ 
tion on the crystal and electronic structure of MeTel 
(Me=Bi, Au, or Pd) compounds, we carried out our cal¬ 
culations at two levels of approximations: the density- 
functional calculations performed within the generalized 
gradient approximation (GCA) using the functional of 
Perdew, Burke, and Ernzerhof (PBE)^° without and with 
the semiempirical dispersion correction. Hence, we com¬ 
pare the results of the PBE (semilocal) and PBE-D2 
(dispersion-corrected) calculations to each other and to 
the experimental data. The total energy calculations per¬ 
formed for compressed and dilated systems are used to 
derive equation of state for metal tellurohalides under 
consideration. Our results show that inclusion of van der 
Waals interactions is necessary for an adequate physical 
description of these compounds. In particular, predicting 
the anisotropy of BiTel via the PBE calculations yields 
an unrealistic result whereas the results of the PBE-D2 
calculations agree with the experimental compressibility 
data. 

Although oxidation state of Bi, Au, and Pd atoms is 
the same (+3) in the BiTel, AuTel, and PdTel com¬ 
pounds, the crystal structures of the latter are quite 
different as shown in Pigs. I(a)-(c). BiTel, AuTel, and 
PdTel crystallize in trigonal, monoclinic, and tetragonal 
structures with space groups P3ml (No. 156), P2i/c 
(No. 14), and PA^/mmc (No. 131), respectively. The 
crystal structure of BiTel could be characterized by the 
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FIG. 1. (Color online) The crystal structures of BiTel (a), AuTel (b), and PdTel (c) compounds. A color contour plot of the 
computed electronic charge density is also shown for each compound. 


hexagonal lattice parameters a and c, and two internal 
parameters since Bi, Te, and I atoms occupy the la, 
Ic, and lb positions with fractional coordinates (0,0,0), 
(2/3,l/3,ZTe), and (1/3,2/3, ^i), respectively, cf. Ref. 15. 
As shown in Figure 1(a), adjacent trigonal layers formed 
by Bi, Te, and I atoms stack along the c-axis of the hexag¬ 
onal lattice in BiTel. The crystal structure of AuTel could 
be characterized by the monoclinic lattice parameters a, 
b, c, and /3, and nine internal parameters since Au, Te, 
and I atoms occupy the 4e positions with fractional co¬ 
ordinates (xau,2/Au,zau), (xTe,PTe,ZTe), (xi,pi,zi), respec¬ 
tively, cf. Ref. 18. As shown in Fig. 1(b), the atoms 
form corrugated layers parallel to the fee-plane, which 
stack along the a-axis in AuTel. The crystal structure 
of PdTel could be characterized by the tetragonal lat¬ 
tice parameters a and c, and three internal parameters 
since Pd, Te, and I atoms occupy the 4m, 41, and 4j 
positions with fractional coordinates (0,ypd,0), (xTe,0,0), 
and (xi, 1/2,0), respectively, cf. Ref. 19. As shown in 
Fig. 1(c), ladder chains made of Te 2 Pdl 2 units running 
along the a- and fe-axes stack along the c-axis in PdTel. 
Note that the Te and I atoms are mixed in the layers or 
ladder chains of AuTel or PdTel, respectively, whereas 
the layers of BiTel are unary, consisting of Bi, Te or I 
atoms. Accordingly, the space group of the BiTel crys¬ 
tal is noncentrosymmetric while the crystal structures 
of AuTel and PdTel are centrosymmetric. This puts 
BiTel in a unique position in terms of spintronic ma¬ 
terial properties,^ which lacks inversion symmetry in the 


TABLE I. Dispersion coefficients Ce (in J nm® mol“^) and 
van der Waals radii Ro (in A) employed in the present PBE- 
D2 calculations. These values are taken from the references 
given in the rightmost column. 


Element 

Ce 

Ro 

Reference 

Pd 

24.67 

1.639 

16 

Te 

31.74 

1.892 

16 

I 

31.50 

1.892 

16 

Au 

40.62 

1.772 

21 

Bi 

63.55 

1.900 

22 


crystal structure. 

The rest of the paper is organized as follows: Sec. II 
is devoted to the method of calculation, giving also a 
summary of the computational details. We focus on 
the crystal structure, equation of state, and electronic 
band structure in Sec. Ill where we discuss the results of 
our semilocal (PBE) or dispersion-corrected (PBE-D2) 
density-functional calculations. A summary of our find¬ 
ings is given in Sec. IV. Lastly, the spin-orbit coupling 
(SOC) effects on the BiTel band structure are described 
in further details in the Appendix. 


II. METHOD 

All calculated properties reported here were obtained 
via semilocal^^ (PBE) or dispersion-corrected^® (PBE- 
D2) density-functional calculations performed by em¬ 
ploying the projector augmented-wave (PAW) method,^® 
as implemented^^ in the Vienna ab initio simulation 
package^® (VASP). The calculations for BiTel and AuTel 
were performed in the noncollinear mode^®’^^ of VASP 
in order to take spin-orbit coupling into account. The 
4d and 5s, 5s and 5p, 5s and 5p, 5d and 6s, and 6s 
and 6p states are treated as valence states for palla¬ 
dium, tellurium, iodine, gold, and bismuth, respectively. 
Plane wave basis sets were used to represent the elec¬ 
tronic states, which were determined by imposing a ki¬ 
netic energy cutoff of 325 eV. The long-range disper¬ 
sion corrections^® for periodic systems were treated as 
described in Ref. 28. The values of the dispersion co¬ 
efficient Cq and van der Waals radius Rq used in this 
study are given in Table I. The global scaling factor sg 
was set to 0.75, which is the adequate value for the PBE 
functional.^® 

We first carried out full optimization of the crystal 
structures where concurrent relaxations of the unit cell 
volume and shape as well as the ionic positions were per¬ 
formed with no symmetry constraints, until the maxi¬ 
mum value of residual forces on atoms was reduced to 
be smaller than 0.01 eV /A. Convergence criterion for the 
electronic self-consistency was set up to 10“® eV. In these 
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optimizations, we used the primitive unit cells of BiTel, 
AuTel, and PdTel, whose Brillouin zones were sampled 
by 20 X 20 X 16, 8 x 8 x 8, and 21 x 21 x 29 k-point 
meshes, respectively, which were generated according to 
Monkhorst-Pack scheme.^® Using the optimized crystal 
structures, we then carried out band-structure and (pro¬ 
jected) density-of-states calculations. Besides, we per¬ 
formed geometry optimizations for the elemental solids 
of bismuth, gold, palladium, tellurium, and iodine, and 
employed the respective equilibrium energies per atom 
Emb (Me=Bi, Au, Pd), E^e and Ei in the computation 
of the formation energy AHf. It should be reminded that 
the form of the electronic Hamiltonian used in dispersion- 
corrected (PBE-D2) calculations is the same as in the cal¬ 
culations employing the PBE functional alone. In other 
words, the effect of dispersion correction on the electronic 
structure is indirectly through modihcation of the crystal 
structure since the van der Waals interactions are treated 
as semiempirical force fields in the DFT-D2 approach.^® 
Secondly, we carried out constant-volume optimization 
of the crystal structures where the unit cell shape and the 
ionic positions were allowed to relax. Hence, we obtained 
the energy E per formula unit as a function of the vol¬ 
ume V per formula unit, which was used to derive equa¬ 
tion of state (EOS) at zero temperature. We found that 
the energy-volume curve is not accurately reproduced 
by a third-order Birch-Murnaghan (BM) fit that is in 
widespread use, which is further discussed in Section HI. 
Thus, we performed forth- and fifth-order BM fits®®’®^ as 
a function of the Eulerian strain / = [(Vq/U)^/® — I]/2 
that is dehned from V and the zero-pressure volume Vg, 
employing 

h 

'^max 

E = Eg+Y, Ck-lf ( 1 ) 

k=2 

with fcinax = 4 and 5, respectively. Here CkS are the 
fitting coefficients, and Eg denotes the equilibrium energy 
(per formula unit). Note that AEIf = Eq — (i?Me + Axe + 
El). The pressure P was computed by using 


P = 


(l + 2/)t 

3Uo 


h — ^ 

'''max 

^ ik + l)C,f. 




( 2 ) 


The isothermal bulk modulus Kg and its pressure deriva¬ 
tives Kq, Kq, and Kq' (all evaluated at Vg) were obtained 
via 
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respectively. 

The volume compressibility Ky was obtained as the in¬ 
verse of the bulk modulus, i.e., Ky = IjKg. The ax¬ 
ial (linear) compressibilities ki = — {dlnl/dP)p^Q, with 
I denoting the lattice constant along one of the crys¬ 
tal axes, were computed by dividing — j by 

iw)v=Vo’ '^bere the former [latter] was obtained via 
cubic spline interpolation of the lattice parameter I as 
a function of the volume V [via Eq. (2)]. The relia¬ 
bility of this procedure was tested by checking whether 
Ky = Ka + Kb + Kc holds Or not. The latter equality was 
satisfied in the cases of BiTel and PdTel, but not in the 
case of AuTel. As discussed in Section HI, the variation 
of the lattice parameter c with the pressure is nonmono¬ 
tonic in the case of AuTel, which results in a substantial 
error in estimating derivatives via spline interpolation. 
The linear compressibilities of AuTel are therefore not 
reported. 


III. RESULTS AND DISCUSSION 
A. Crystal Structure 


The crystal structure optimizations result in the ex¬ 
perimentally determined ground-state structures, i.e., a 
noncentrosymmetric (centrosymmetric) crystal for BiTel 
(AuTel and PdTel) that are shown in Fig. l(a)-(c) where 
color contour plots of computed electronic charge density 
are also displayed. Inspection of the charge-density plot 
in Fig. 1(a) reveals that the most electron-rich (electron- 
poor) regions in BiTel are around I (Bi) atoms. This 
reflects the fact that bismuth is less electronegative than 
both iodine and tellurium.®^ It is also interesting to 
note that the electronegativity difference Xxe ~ XBi is 
positive and comparatively small whereas xi ~ XBi is 
also positive but comparatively large.^"^'^^ Accordingly, 
the electronic charge distribution around Bi—Te (Bi—I) 
bonds in Fig. 1(a) is rather of a covalent (ionic) charac¬ 
ter, which introduces an asymmetry between Bi—Te and 
Bi—I bonds. On the other hand, we have yxe — XMe < 
0 < XI - XMe < Xi - XTe for both Me=Au and Pd, i.e., 
the differences XTe — XMe and xi ~ XMe are both rela¬ 
tively small, which are of the opposite sign.^^iSS agree¬ 
ment with the latter, the electronic charge distribution 
around not only Me—Te bonds and but also Me—I bonds 
in Fig. l(b)-(c) look more like those of covalent bonds. 
Thus, the tellurium and iodine atoms prefer to coordinate 
with the metal atoms almost equally, and therefore tend 
to mix (as opposed to form unary layers) in the AuTel 
and PdTel compounds, rendering the crystal structures 
of the latter centrosymmetric. 

A comparison of the results of the crystal structure 
optimizations (PBE and PBE-D2) to the experimental 
data is given in Tables H and SI (Ref. 34). It is seen 
in Table H that the improvement due to the dispersion 
correction is mostly on the lattice parameters (a, b, c) 
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TABLE II. Calculated (PBE and PBE-D2) and measured lat¬ 
tice parameters of the MeTel crystals. The experimental val¬ 
ues are taken from Ref. 15, Ref. 18, and Ref. 19 for BiTel, 
AuTel, and PdTel, respectively. 




PBE 

PBE-D2 

Exptl. 

BiTel 


a (A) 

4.4371 

4.2843 

4.3392 


c (A) 

7.433 

7.021 

6.854 


ZTe 

0.7692 

0.7479 

0.6928 


Zl 

0.2828 

0.3115 

0.2510 

AuTel 


a (A) 

8.0057 

7.2579 

7.3130 


b(A) 

7.8918 

7.3654 

7.6242 


C (A) 

7.4208 

7.3483 

7.2550 


P 

104.81 

103.95 

106.26 


XA^ 

0.4589 

0.4746 

0.4654 


VAu 

0.1333 

0.1370 

0.1395 


ZAu 

0.2308 

0.2405 

0.2370 


XTe 

0.6481 

0.6868 

0.6720 


VTe 

0.1433 

0.1331 

0.1301 


ZTe 

0.9821 

0.9908 

0.9910 


Xl 

0.1954 

0.1626 

0.1758 


yi 

0.0659 

0.1059 

0.0857 


Zl 

0.3895 

0.3799 

0.3892 

PdTel 


a (A) 

8.095 

7.806 

7.821 


c (A) 

5.722 

5.652 

5.659 


ypd 

0.2500 

0.2495 

0.2525 


XTe 

0.2086 

0.2262 

0.2164 


Xl 

0.2321 

0.2462 

0.2435 


whereas the PBE and PBE-D2 calculations yield errors 
of similar magnitude in the prediction of the internal pa¬ 
rameters (x, y, z). It is thus notable that the prediction 
of unit cell volume and shape is substantially improved: 
Figure 2(a) displays a plot of the calculated (PBE and 
PBE-D2) versus experimental values for the equilibrium 
volume per formula unit. Although it is well known the 
unit cell volume is overestimated within the GGA, the er¬ 
ror in the PBE-optimized volume is clearly greater than 
expected, which is 13.4, 16.7, and 8.3 % for BiTel, AuTel, 
and PdTel, respectively. As evident from the trend of the 
empty symbols in Fig. 2(a) this anomalous overestima¬ 
tion is widespread. Since the lattice parameters (a, b, c) 
obtained in the PBE-D2 optimizations are smaller than 
those in the PBE-D2 optimizations, cf. Table II, the pre¬ 
diction of equilibrium volume is systematically improved 
thanks to the dispersion correction, as evident from the 
trend of the filled symbols in Fig. 2(a). Moreover, a com¬ 
parison of the empty (PBE) and filled (PBE-D2) symbols 
of the same shape to each other in Fig. 2(b) shows that 
the errors in the prediction of b/a and c/a ratios are, at 
the same time, significantly reduced thanks to the dis¬ 
persion correction. 

As for the prediction of the internal parameters of the 
BiTel crystal, it is interesting to compare our results to 
the results of Bahramy et al.^ (Sklyadneva et al.^) who 
predicted zxe = 0.7482 and zi = 0.3076 (^Te = 0.758 
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FIG. 2. (Color online) The plots of calculated versus experi¬ 
mental equilibrium volume (a), and error in b/a or c/a ratios 
(b). 
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and zi = 0.299). Note that our PBE-D2 and PBE values 
for ZTe and zi are in close agreement with those given 
by Bahramy et al.^ and Sklyadneva et al.,^ respectively 
-albeit in both Refs. 5 and 8 the PBE functional was em¬ 
ployed with no additional terms for the dispersion cor¬ 
rection. This is so because the lattice parameters a and 
c were fixed to their experimental values (that are close 
to our PBE-D2 values) in Ref. 5, which was not done in 
Ref. 8. 

The predictions of Ref. 5 and 8 show that the computed 
values of the equilibrium Bi—Te and Bi—I distances cor¬ 
respond to the Bi—I and Bi—Te distances in the exper¬ 
imentally determined structure,^® respectively, which is 
supported by our PBE-D2 values, cf. Table SI (Ref. 34). 
Note that in Ref. 15 the bond lengths given in Table 3 are 
not consistent with the internal parameters given Table 2. 
It is therefore problematic to perform a direct compari¬ 
son of experimental and optimized values of the internal 
lattice parameters zxe and zj. Thereupon, we consid¬ 
ered a second “phase” for BiTel, which was generated 
by exchanging the positions of Te and I atoms so that 
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the Te and I atoms occupy the 16 and Ic positions with 
fractional coordinates (l/3,2/3,2;Te), and (2/3,1/3, 2 : 1 ), re¬ 
spectively. The crystal structure optimization performed 
for this second phase yielded zxe = 0.2520 and zi = 
0.6884, and the total energy the same as that of the BiTel 
phase described in Sec. I and Table II. In line with the 
latter, comparison performed with the aid of the COMP- 
STRU program^® proved that this second BiTel phase 
(Te atoms residing at 16 positions with ZTe = 0.2520 and 
I atoms residing at 16 positions with zi = 0.6884) is in¬ 
deed identical to the first one (Te atoms residing at Ic 
positions with ZTe = 0.7479 and I atoms residing at 16 
positions with zi = 0.3115). It is interesting to point out 
that excellent agreement between experimental data and 
our computed values (zxe = 0.2520 and Zj = 0.6884) is 
obtained once the positions of Te and I atoms are ex¬ 
changed in the experimental^® crystal structure of BiTel. 
As long as the x-ray diffraction fails to distinguish the Te 
and I layers in BiTel,® we anticipate that a full-fledged 
experimental characterization would render the values of 
ZTe and zi in agreement with our optimized (PBE-D2) 
values given in Table II. 


B. Equation of State 

A plot of the energy difference AE = E — (Eue + 
Etc + El) with respect to the volume V is given in 
Figs. 3(a), (b), and (c) for BiTel, AuTel, and PdTel, 
respectively. The calculated values are represented by 
the blue (PBE) and red (PBE-D2) circles. The solid¬ 
line curves connecting the symbols show the forth- and 
fifth-order BM fits, whose equation is given by AE = 
/k, cf. Eq. (I). Although we first 
performed the BM fits with fcmax = 3, we found it neces¬ 
sary to increase fcmax to 4 for AuTel and PdTel and to 5 
for BiTel in order to ensure a satisfactory level of accu¬ 
racy in the fitting procedure, as mentioned in Section II. 
It is discernible in Fig. SI (Ref. 34) that a good fit is not 
obtained when fcmax is reduced to 4 for BiTel, and 3 for 
AuTel and PdTel. 

The values of Ck parameters obtained via fitting as 
well as the formation energies AElf are given in Ta¬ 
ble III. Note that while AElj is significantly (AuTel) 
or slightly (PdTel) reduced due to addition of the dis¬ 
persion terms to the PBE functional, both PBE and 
dispersion-corrected (PBE-D2) calculations yield a sim¬ 
ilar degree of overbinding for BiTel since the experi¬ 
mental value®® of the enthalpy of formation for BiTel 
is Ai7®(298 K) = 0.856 eV per formula unit. On the 
other hand, the interlayer binding energy of BiTel was 
computed to be 24 and 332 meV per formula unit in our 
PBE and PBE-D2 calculations, respectively. Thus, the 
dispersion-corrected (PBE-D2) calculations yield a sub¬ 
stantially stronger interlayer binding for BiTel, compared 
to the PBE calculations. It is also notable that the PBE- 
calculated binding curves of BiTel [Fig. 3(a)] and AuTel 
[Fig. 3(b)] are anomalously flat in the vicinity of equilib¬ 


rium volume Vo (and especially for volumes larger than 
Vo), which is not the case for PdTel [Fig. 3(c)]. This flat¬ 
ness of the binding curves could be attributed to lack¬ 
ing van der Waals interactions because the dispersion- 
corrected (PBE-D2) binding curves have a significantly 
increased curvature. 

The pressure-volume curves are shown in Figs. 3(d), 
(e), and (f) for BiTel, AuTel, and PdTel, respectively. 
The graphs given in Figs. 3(g)-(m) shows the variation 
of the lattice parameters (a, 6, c) with the pressure. The 
experimental compressibility data® is also included in 
Figs. 3(d), (g), and (h) for comparison. The experimen¬ 
tal data employed in Fig. 5 of Ref. 14 is not presented 
in Fig. 3(d) since is seems to be shifted in comparison 
to the data of Ref. 3, cf. Fig. S2 (Ref. 34), which also 
does not agree with the equilibrium volume reported in 
Ref. 15. It is noticeable in Fig. 3(d) that the dispersion- 
corrected V (P) curve is in much better agreement with 
the experimental data (the empty circles) whereas the 
PBE-calculated curve has a substantially larger slope in 
the low-pressure region. The latter appears to be the 
case with AuTel too, cf. Fig. 3(e). In contrast, the PBE- 
calculated and dispersion-corrected (PBE-D2) curves for 
PdTel have comparable slopes even in the low-pressure 
region. Note that the steeper decay of V{P) curve in the 
low-pressure region could be traced back to the variation 
of lattice parameters c and a for BiTel and AnTel, re¬ 
spectively, cf. Figs. 3(h) and (i). Thereupon, the error 
in the PBE calculations is clearly not restricted to the 
overestimation of the volume for quasi-layered systems, 
e.g.. Figs. 1(a) and (b). One should, on the other hand, 
also notice that the agreement between the PBE-D2 and 
experimental curves in Fig. 3(d) is rather coincidental 
since the lattice parameter c (a) is overestimated (under¬ 
estimated) in the PBE-D2 calculations, regardless the 
value of P, cf. Figs. 3(g) and (h). Furthermore, the cal¬ 
culated and experimental values in Fig. 3 refer to zero 
and room temperature, respectively. Therefore, a com¬ 
parison of the relative volume {V/Vo) versus normalized 
pressure (P/Kq) curves is given in Fig. 4 for BiTel, fol¬ 
lowing a proposal of Ref. 37 [where using a simple scaling 
of variables was demonstrated to be practical in com¬ 
paring the calculated P(V) equation to the experimental 
equation of state]. It is seen in Fig. 4 that the experimen¬ 
tal data points are bracketed by the PBE-calculated and 
dispersion-corrected (PBE-D2) curves, which are seem¬ 
ingly in better agreement with the PBE-D2 curve than 
the PBE curve. 

Turning back to Fig. 3, it is interesting to point out 
that the variation of the lattice parameter c of AuTel 
with the pressure is slightly increasing (nonmonotonic) 
according to the results of the PBE-D2 (PBE) calcu¬ 
lations, as shown in Fig. 3(k). This implies that the 
lattice of AuTel would expand along the c-axis under 
compression. In other words, the linear compressibility 
Kc of AuTel is negative according to our calculations, cf. 
Fig. 3(k). It should be remarked that this prediction calls 
for experimental verification inasmuch as the compress- 
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FIG. 3. (Color online) The plot of the energy difference AE — E — {Emb + Etb + Ei), cf. Eq. (1), versus the volume for BiTel 
(a), AuTel (b), and PdTel (c). Here fcmax is set to 4 for BiTel, and 3 for AuTel and PdTel. The respective pressure-volume 
curves are plotted in (d), (e), and (f). The variation of the lattice parameters (a, h, c) with the pressure is shown in (g)-(m). 
The PBE-calculated and dispersion-corrected (PBE-D2) curves are in blue and red, respectively. The filled symbols mark the 
experimental values of the equilibrium volumes [(d)-(f)] and lattice parameters [(g)-(m)] measured by Shevelkov et al. (Ref. 15), 
Fenner et al. (Ref. 18), and Seo et al. (Ref. 19) for BiTel, AuTel, and PdTel, respectively. The empty circles in (d), (g), and 
(h) represent the experimental compressibility data of Xi et al. (Ref. 3) for BiTel. 


TABLE III. The values of the formation energy AHf, the equilibrium volume Vo, and the fitting coefficients Ck in Eq. (1) for 
MeTel (Me=Bi,Au,Pd) compounds. 




h 

•^max 

AHf (eV) 

Vo(k^) 

Cl (eV) 

Ga (eV) 

Co (eV) 

G 4 (eV) 

BiTel 

PBE 

5 

-0.958 

127.039 

15.357 

278.613 

1091.512 

-8932.044 


PBE-D2 

5 

-0.978 

111.626 

88.806 

228.317 

-4645.020 

30756.053 

AuTel 

PBE 

4 

-0.404 

113.317 

13.558 

172.994 

-444.287 



PBE-D2 

4 

-0.337 

95.308 

58.526 

244.128 

-1103.027 


PdTel 

PBE 

4 

-0.973 

93.742 

61.728 

-11.560 

506.408 



PBE-D2 

4 

-0.961 

86.101 

90.574 

103.902 

-695.661 
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FIG. 4. (Color online) The plot of the relative volume V/Vb 
versus the normalized pressure PjKo for BiTel. The PBE- 
calculated and dispersion-corrected curves are in blue and 
red, respectively. The black (solid) circles represent the ex¬ 
perimental data of Xi et al. (Ref. 3). The black curve was 
obtained by performing a third-order BM fit to the latter. 


ibility of AuTel has not been investigated before, to our 
knowledge. Note that the usual behavior under compres¬ 
sion is that the individual lattice parameters decrease so 
that Ki > 0, which is the case in Figs. 3(g)-(j) and (1)- 
(m). The behavior of negative linear compressibility has 
nevertheless been observed in a number of systems^®^"^^ 
recently. 

The bulk modulus Kq and its pressure derivatives Kq, 
Kq, and Kq' computed via Eq. (3) with Ck coefficients 
listed in Table III are given in Table IV where the ex¬ 
perimental values for BiTel are also included. The lat¬ 
ter were obtained by performing a third-order BM fit, 
cf. the black solid curve in Fig. 4, to the experimen¬ 
tal compressibility data.® It is seen that the PBE cal¬ 
culation yields a substantially underestimated (overesti¬ 
mated) value for Kq {Kq) of BiTel. Recalling that (i) the 
PBE-calculated binding energy curve anomalously flat in 
the vicinity of equilibrium volume, cf. Fig. 3(a), and (ii) 
the PBE-calculated equation of state has an excessively 
steep slope in the low-pressure region, cf. Fig. 3(d), it is 


TABLE IV. The bulk modulus Kq (in GPa) and its pressure 
derivatives Kq, Kq (in GPa“^), and Kq' (in GPa“^). 




^max 

Kq 

K'q 

K 

K' 

BiTel 

PBE 

5 

4.3 

22.1 

-59.6 

658.8 


PBE-D2 

5 

28.3 

6.6 

-2.9 

2.5 


Experimental 

3 

20.5 

7.6 



AuTel 

PBE 

4 

4.3 

16.8 

-52.4 



PBE-D2 

4 

21.9 

8.2 

-2.3 


PdTel 

PBE 

4 

23.4 

3.8 

0.3 



PBE-D2 

4 

37.5 

5.1 

-0.4 



TABLE V. The volume (k„) and axial {na and Kc) compress¬ 
ibilities (in GPa“^), and their ratios. 




Kv 

Ka 

f^c 

KajKy 

KcjKy 

KajKc 

BiTel 


PBE 

0.232 

0.010 

0.212 

0.043 

0.914 

0.047 


PBE-D2 

0.035 

0.011 

0.013 

0.314 

0.371 

0.846 

PdTel 


PBE 

0.043 

0.017 

0.008 

0.395 

0.186 

2.125 


PBE-D2 

0.027 

0.010 

0.007 

0.370 

0.259 

1.429 


no surprise that the bulk modulus is substantially under¬ 
estimated in the PBE calculations. The PBE-calculated 
first pressure derivative K'q is consequently greatly over¬ 
estimated, balancing this underestimation of Kq. On 
the contrary, the dispersion-corrected (PBE-D2) calcu¬ 
lations, albeit in much better agreement with the experi¬ 
mental values, result in overestimation (underestimation) 
for Kq (Kq). Hence, the experimental values of Kq and 
K'q are bracketed by the PBE-calculated and dispersion- 
corrected (PBE-D2) values. 

The volume (bulk) and axial (linear) compressibili¬ 
ties and their ratios are given in Table V for BiTel and 
PdTel. Note that Kb = Ka for these systems. For BiTel 
(PdTel), the PBE and PBE-D2 values for the a-axis (c- 
axis) compressibility are close to each other. In con¬ 
trast, K^^^ < for BiTel, and 

for PdTel. Thus, the c-axis compressibility of BiTel is 
greatly overestimated in the PBE calculations (which re¬ 
flects lacking van der Waals binding). We have further¬ 
more KajK q <C 1 (PBE) while KajK c ~ 1 (PBE-D2), and 
i^al^y < 1/3 and Kc/kv > 1/3 (PBE) while KajKy ^ 1/3 
and KcIkv ^1/3 (PBE-D2) for BiTel. Note that the ra¬ 
tio KajKc should be unity, and both Kafn^ and KqIkv 
need to be equal to 1/3 for isotropic materials. For lay¬ 
ered materials, on the other hand, one would expect to 
have the ratio Kaj Ky {kq/K y) significantly smaller (larger) 
than 1/3 and the ratio KajKc much smaller than unity 
(e.g., Ka/Ky = 0.028, Kc/Ky = 0.943, KajKc = 0.03 for 
graphite,a most known layered material). Thus, the 
PBE calculations predict BiTel to be an extremely lay¬ 
ered material whereas the PBE-D2 calculations indicate 
a layered, but substantially less anisotropic material. In 
other words, predicting the anisotropy of BiTel via the 
PBE calculations yields an unrealistic result, given that 
the PBE-D2 calculations are in much better agreement 
with the experimental compressibility data. It appears 
that the inclusion of van der Waals attraction, which was 
motivated by the quasi-layered structure of BiTel, leads 
to a much less anisotropic crystal structure. Note that 
this issue does not raise in the case of PdTel where van 
der Waals binding plays a much less prominent role and 
therefore the PBE and PBE-D2 calculations qualitatively 
yield similar results, cf. Table V. 













C. Electronic Structure 


Figure 5(a)-(c) displays the calculated band structure 
(left panels in each subplot) and projected density of 
states (right panels in each subplot) of the MeTel com¬ 
pounds under consideration, obtained by using the op¬ 
timized (PBE-D2) lattice parameters given in Table 11. 
BiTel and AuTel appear to be narrow band gap semicon¬ 
ductors whereas PdTel is predicted to be a metal with 
relatively low density of states at the Fermi level Ep. In 
Fig. 5(a), Rashba-type splitting near the Brillouin zone 
point A is visible, which is in agreement with the former 
studies.^’® This is important because the emergence of 
the bulk Rashba-type splittings in the calculated band 
structure depends on the reliable prediction of the BiTel 
lattice parameters, cf. Fig. 4 of Ref. 13. We present a 
comparison of the BiTel band structures with and without 
SOC in the Appendix, which confirms that splitting near 
the point A in Fig. 5(a) manifests SOC-induced modifica¬ 
tion of the band edges. It is further shown that the gra¬ 
dient of the crystal potential VR(r) leads predominantly 
to an asymmetric electric field along the c-axis that facil¬ 
itates the spin-orbit coupling. As long as the latter gives 
rise to an effective magnetic field"^^ Hgg, a simple under¬ 
standing is that Rashba-type splitting in Fig. 5(a) orig¬ 
inates from the interaction of the crystal electrons with 
HeS that is proportional to VR (r) x p, where p denotes 
the momentum operator. It is to be noted that a compre¬ 
hensive analysis has already been given by Bahramy et 
al. in Ref. 5, revealing the origin of (giant) Rashba-type 
splitting in BiTel. 

We see in Fig. 5(a) that the lower valence band con¬ 
sists of the Bi, Te, and I s orbitals whereas the upper 
valence band as well as lower conduction band are de¬ 
rived from the Bi, Te, and I p orbitals. Fig. 5(b) shows 
that the bands of AuTel are relatively flat (nondisper- 
sive), which are derived primarily from Te and I s orbitals 
(lower valence band), and Te and I s and Au d orbitals 
(upper valence band as well as lower conduction band). 
Au d states dominate the bottom of the upper valence 
band. Fig. 5(c) shows that the bands of PdTel around 
the Fermi level are quite dispersive along not only per¬ 
pendicular F—Z direction but also parallel F-M and X—M 
directions. Thus our calculations indicate that opening a 
band gap owing to a possible c-axis-doubling distortion^® 
would probably not render PdTel a semiconductor. We 
see in Fig. 5(c) that the deep-lying bands consists mainly 
of the Te s and I s orbitals while the bands in the vicinity 
of the Fermi level are derived from the Te p, I p and Pd 
d orbitals. 

Table VI gives some characteristic band structure pa¬ 
rameters for BiTel and AuTel semiconductors, for which 
experimental data is available: the band gap Eg and the 
Rashba energy Ep and momentum offset kp for the con¬ 
duction band minimum. As for the BiTel band gap, 
it is interesting to point out that our PBE-D2 value is 
close the PBE value (0.242 eV) of Rusinov et al}'^ who 
employed the experimental^® lattice parameters a and c 


TABLE VI. The calculated (PBE and PBE-D2) and exper¬ 
imental band structure parameters of the BiTel and AuTel 
semiconductors: the band gap Eg (in eV) and the Rashba 
energy Ep (in eV) and momentum offset kp (A“^) for the 
conduction band minimum. 




PBE 

PBE-D2 

Experimental 

BiTel 

E, 

0.432 

0.207 

0.38“, 0.36\ 0.26“ 


Ep 

0.081 

0.124 

0.1“ 


kp 

0.0529 

0.0514 

0.052“ 

AuTel 

Eg (indirect) 

0.938 

0.623 



Eg (direct) 

0.938 

0.803 

O.O'' 

“Ref. 1 

*'Ref. 44 

“Ref. 45 

‘'Ref. 

17 


and PBE-optimized® internal parameters ztb and zi in 
their calculations. Similarly, our PBE-D2 values for the 
Rashba energy and momentum offset are in close agree¬ 
ment with the respective PBE values reported in Ref. 10 
(where Ep =0.122 eV and kp = 0.050 A“^). 

It is clear from the entries of Table VI that the calcu¬ 
lations with the PBE functional yield band gaps that are 
greater than the respective experimental values, which is 
against the expected trend inasmuch as the band gap is 
known to be underestimated within the GGA. Contrary 
to this, the calculated band gaps corresponding to the 
PBE-D2 optimized crystal structures are smaller than 
the respective experimental values, in line with the ex¬ 
pected trend. Hence, the experimental band gaps are 
bracketedhy the respective PBE and PBE-D2 values, i.e.. 


E 


PBE-D2 


< Eg < Eg 


iPBE 


Similarly, the experimental values of Ep and kp are also 
bracketed by the respective PBE and PBE-D2 values (cf. 
Table VI) since 


E^^ < Ep < £;Pbe-D2 


fcPBE-D2 ^ kp <k 


R 

PBE 
R ■ 


We accordingly anticipate that comparative usage of the 
PBE and PBE-D2 calculations would be a useful strat¬ 
egy for predicting the band structure of layered polar 
semiconductors, which eliminates the need to employ the 
experimental lattice parameters. 

As for the Rashba parameter ap = 2Ep/kp of BiTel, 
which represents the strength of Rashba splitting, the 
PBE-calculated and dispersion-corrected (PBE-D2) val¬ 
ues are ap = 3.06 and 4.82 eVA, respectively, bracketing 
the experimental value^ of ap = 3.85 eVA. It is shown 
in the Appendix that ending up with a larger Rashba 
parameter within the dispersion-corrected (PBE-D2) de¬ 
scription is due to the enhancement of the asymmetric 
electric field (—VR) along the c-axis [Fig. A2(b)]. 

Given that the experimental values of kp, Ep, and 
ap are bracketed by their calculated (PBE and PBE- 
D2) values, one might attempt to interpolate their PBE- 
calculated and dispersion-corrected (PBE-D2) values in 
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T X M r Z R A Z 0 10 22 
DOS 


FIG. 5. (Color online) The calculated band structure (left panels in each subplots) and projected density of states (right panels 
in each subplots) of BiTel (a), AuTel (b), and PdTel (c). The crystal structures optimized via dispersion-corrected (PBE-D2) 
calculations are assumed, cf. Table II. 



0.20 0.25 0.30 0.35 0.40 0.45 


(eV) 

FIG. 6. (Color online) The momentum offset kn (a), Rashba 
energy Eh (b), and Rashba parameter or (c) for the conduc¬ 
tion band minimum of BiTel versus the band gap Eg. 


order to obtain better estimates. The band gap Eg can 
be used as an interpolation variable for this purpose since 
the values of kfi and Er are in almost linear correlation 
with the band gap values, as shown in Figs. 6(a) and (b), 
respectively. Note that the error due to the underestima¬ 
tion (PBE-D2) or overestimation (PBE) of Eg could then 
be corrected by using the experimental value of the band 


gap (i.e., Eg = 0.36 eV) in the interpolation. The dashed 
lines in Figs. 6(a) and (b) passes though the calculated 
points, whose equations are given by 


kR = 0.05002 -bO.OObeTEg, 

Efl = 0.16356-O.igillEg, (4) 


respectively. Combining the latter two equations, ur 
could be obtained as a function of Eg, which is shown 
by the dashed line in Fig. 6(c). Substituting the experi¬ 
mental value of Eg in Eq. (4) yields the interpolated val¬ 
ues of kR = 0.0524 A“^, Er =0.095 eV, and aR = 3.62 
eVA [which are marked by the black circles in Figs. 6(a), 
(b) and (c), respectively] in close agreement with the re¬ 
spective experimental values^ [marked by the green cir¬ 
cles in Figs. 6(a)-(c)], cf. Table VI. 

Finally, forasmuch as the pressure-induced closure of 
the BiTel band gap has been explored in a number of re¬ 
cent studies,we estimated the critical pressure 
Pc at which the band gap disappears (the latter being a 
signature of the topological phase transition). We obtain 
Pc = 3.73 GPa using our dispersion-corrected (PBE-D2) 
equation of state, cf. Section IIIB, together with a series 
of band structure calculations (not shown), which con¬ 
firms the refined value {Pc = 3.5 GPa) of Bahramy et 
al. (Refs. 2 and 46). Our investigation of the variation of 
the BiTel band gap and Rashba parameters with pressure 
will be reported in detail in a separate publication. 
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IV. CONCLUSION 

In summary, our comparative investigation of the re¬ 
sults of semilocal (PBE) and dispersion-corrected (PBE- 
D2) calculations reveals the effect of van der Waals at¬ 
tractions on the crystal and electronic structure of metal 
tellurohalides BiTel, AuTel, and PdTel. We find that the 
prediction of crystal structure is systematically improved 
thanks to the inclusion of van der Waals dispersion. It 
is shown for the compounds with a quasi-layered crys¬ 
tal structure, viz. BiTel and AuTel, that (i) the PBE- 
calculated energy versus volume curve is anomalously flat 
in the vicinity of equilibrium volume, and especially for 
volumes larger than the equilibrium volume, and (ii) the 
PBE-calculated equation of state has an excessively steep 
slope in the low-pressure region, which are also fixed 
in the dispersion-corrected (PBE-D2) calculations. Our 
analysis based on the computation of the volume and ax¬ 
ial compressibilities shows that predicting the anisotropy 
of BiTel via the semilocal calculations yields an unreal¬ 
istic result whereas the results of the PBE-D2 calcula¬ 
tions agree with the experimental compressibility data. 
Our calculations render that BiTel and AuTel are nar¬ 
row band gap semiconductors with Rashba-type spin¬ 
splitting at the band edges and with an indirect band 
gap, respectively. PdTel, on the other hand, is predicted 
to be a metal with relatively low density of states at the 
Fermi level. It is notable that the band gaps computed 
by using the optimized crystal structures from the PBE 
(PBE-D2) calculations are found to be greater (smaller) 
than the respective experimental values, which is against 
(in line with) the expected trend inasmuch as the band 
gap is known to be underestimated within the GGA. We 
also find that the Rashba parameter, Rashba energy and 
momentum offset of BiTel are bracketed by the respective 
values obtained via the semilocal (PBE) and dispersion- 
corrected (PBE-D2) calculations. Specifically, a larger 
value for the Rashba parameter is obtained in the PBE- 
D2 calculations, which could be attributed to the reduc¬ 
tion of the band gap caused by modification of the crystal 
structure owing to the inclusion of van der Waals disper¬ 
sion terms. Excellent agreement with the experimental 
data for the Rashba parameter, Rashba energy and mo¬ 
mentum offset of BiTel is obtained via interpolation of 
the calculated (PBE and PBE-D2) values. 
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FIG. Al. (Color online) The band structure of BiTel in the 
vicinity of the band edges, calculated with and without SOC 
and dispersion correction. The lower conduction (upper va¬ 
lence) band states labeled as CBM and CBM-l-1 (VBM and 
VBM—1) are marked by the solid circles. 


Appendix: Spin-Orbit-Induced Modification of the 
BiTel Band Edges 

A comparison of the calculated band structures 
with and without SOC has already been presented by 
Bahramy et al. in Ref. 5 where the PBE-optimized val¬ 
ues for zxe and zi were employed together with the ex¬ 
perimental values for a and c. Here we present a sim¬ 
ilar comparison in Fig. Al where no experimental data 
are used. The lower conduction (upper valence) band 
states are marked by the solid circles labeled as “CBM” 
and “CBM-bl” (“VBM” and “VBM-1”) in Fig. Al. Ta¬ 
ble A gives the projected wave function character of these 
states. It is seen in Fig. Al that Rashba-type splitting 
is absent in the band structures calculated without SOC 
(marked as “PBE-D2 w/o SOC” and “PBE w/o SOC”) 
in both cases with and without the dispersion correction. 
When SOC is taken into account, Rashba-type splitting 
occurs on both (conduction and valence) band edges, cf. 
the curves marked as “PBE-D2” and “PBE” in Fig. Al. 
It is thus clear that Rashba-type splitting in Fig. Al [and 
Fig. 5(a)] originates from spin-orbit interaction. Thanks 
to the latter, the band edge states have nonvanishing con¬ 
tributions from px and Py orbitals in comparison to the 
cases without SOC, which could be seen via comparative 
inspection of rows in Table A. Note, for example, that 
when SOC is ignored the character of CBM wave function 
ipc is mostly of Bi-p^ by 72 % (PBE-D2) [69% (PBE)] and 
contribution from px and Py orbitals is negligible. On the 
other hand, the contribution from Bi-p^ is reduced to 35 
% (PBE-D2) [28% (PBE)] at the expense of nonvanish- 
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TABLE A. The projected wave function character of the lower conduction (CBM and CBM+1) and upper valence (VBM 
and VBM—1) band states, cf. Fig. Al. The entry values represent, in %, the angular-momentum-resolved (s, Px, Py, Pz) 
contributions from all (Bi, Te, or I) atoms. The latter are rounded within 1%, and the contributions less than 1% are ignored. 




Bi 




Te 




I 




State 


S 

Px 

Py 

Pz 

s 

Px 

Py 

Pz 

S 

Px 

Py 

Pz 

CBM-bl 

PBE-D2 

2 

13 

10 

22 

2 

14 

15 

4 

1 



14 


PBE 


11 

9 

30 

3 

14 

17 

1 

2 

1 


9 

CBM 

PBE-D2 w/o SOC 




69 

12 




6 



11 


PBE-D2 


12 

16 

28 

5 

14 

10 


2 

2 

3 

6 


PBE 


10 

13 

35 

5 

13 

8 


2 

2 

5 

4 


PBE w/o SOC 




72 

10 




6 



9 

VBM 

PBE-D2 w/o SOC 

16 



1 


3 

9 

44 




24 


PBE-D2 

14 



2 


6 

12 

37 


1 

1 

23 


PBE 

11 



4 


9 

18 

34 


1 

1 

16 


PBE w/o SOC 

13 



5 


8 

17 

33 


2 

4 

16 

VBM-1 

PBE-D2 

13 

4 

2 

1 


7 

10 

35 

1 

6 

3 

14 


PBE 

10 

4 

2 



13 

14 

32 


5 

4 

13 


ing contributions from and Py orbitals when SOC is 
taken into account. Thus, as demonstrated in Fig. A2(a), 
the character of CBM wave function would facilitate for¬ 
mation of a two-dimensional electron gas once the Fermi 
level is above the conduction band minimum. Note that 
the latter is usually the case with the samples used in 
recent experimental studies, e.g., Refs. 1, 12, and 47. It 
is seen in Fig. A2 that the conduction electrons are really 
subject to an asymmetric electric field (—VF) along the 
c-axis [Fig. A2(b)] as mentioned in Section IIIC, which 
are to a signihcant extent confined in the interplanar re¬ 
gions (i) between Te and Bi layers as well as (ii) between 
the Bi and I layers [Fig. A2(a)]. 

Figure A2(b) also shows that the asymmetric electric 
field along the c-axis is substantially enhanced within 
the dispersion-corrected (PBE-D2) description (the red 
curve), compared to the PBE-calculated (blue) curve. 
Note, for comparison, that there is an insignificant 
change in the ab planar average of VF, cf. Fig. S3 
(Ref. 34). Clearly, the enhancement of the asymmetric 
electric field along c-axis leads to stronger SOC. Hence 
the value of Rashba parameter an obtained via the 
dispersion-corrected (PBE-D2) calculation turns out to 
be greater than the respective PBE-calculated value. 

Going back to Fig. Al, it is interesting to point out 
that BiTel turns out to be an indirect semiconductor in 
all calculations presented in Fig. Al, which is in line with 
the former studies.The difference fccBM — fcvBM 
is largest [0.181 A"! for PBE w/o SOC[ when SOC 
and dispersion correction are both ignored, which is re¬ 
duced to 0.044 (0.077 A“^) via inclusion of SOC 

(dispersion correction). Taking SOC and dispersion cor¬ 
rection into account together yields the smallest differ¬ 
ence fccBM — kvBM = 0.010 A“^ (PBE-D2). Similarly 
the band gap Eg is reduced from 1.313 eV (PBE w/o 



z{A) 

FIG. A2. (Color online) The ab planar average of |'!/)c(r)p (a) 
and VF (r) (b) as a function of distance 2 along the c-axis of 
BiTel, obtained via semilocal (PBE) and dispersion-corrected 
(PBE-D2) calculations. The location of the unary Bi, Te, and 
I layers are marked by the vertical dashed lines. 


SOC) to 1.246 eV (PBE-D2 w/o SOC) owing to the dis¬ 
persion correction, which is further reduced to 0.207 eV 
(PBE-D2) via inclusion of SOC. Thus, the reduction in 
the value of Eg is largely due to SOC (rather than the 
dispersion correction). Accordingly, strong (weak) SOC 
implies a smaller (larger) band gap, which corresponds 
to a greater (smaller) value of the Rashba parameter. 
Hence Eg and an are roughly inversely proportional to 
each other, which is also noted in a recent study^° on the 
(BiTeI)m(Bi 2 Te 3 )„ heterostructures. 
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SUPPLEMENTAL MATERIAL 


Table SI lists the calculated and experimental bond 
lengths and bond angles in the MeTel (Me=Bi, An, 
Pd) crystals. 

Figure SI displays a plot of the energy differ¬ 
ence AE versus the volume, and the BM fits with 
fcmax = 4 for BiTel and fcmax = 3 for AuTel and 
PdTel. 

Figure S2 shows a comparison of the calculated and 
experimental pressure-volume curves for BiTel. 

Figure S3 displays the be planar average of the gra¬ 
dient of the BiTel crystal potential as a function of 
distance x along the a-axis. 



S-2 


TABLE SI. The calculated (PBE and PBE-D2) and experimental bond lengths d (in A) and bond angles 6 (in °) of the MeTel 
crystals. The experimental values are taken from Refs. 15, 18, and 19 for BiTel, AuTel, and PdTel, respectively. The labeling 
of bond lengths and angles are therefore the same as in the latter references. 




PBE 

PBE-D2 

Experimental 

BiTel 


dui-Te 

3.08 

3.04 

3.27“ 


dm-i 

3.31 

3.30 

3.04“ 


dxo-Bi-Te 

92.03 

89.55 

83.05*’ 


dxc-Bi-I 

174.45 

174.10 

174.43*’ 


dxo-Bi-I 

91.82 

94.63 

92.79*’ 


dl-Bi-I 

84.06 

80.90 

91.11*’ 

AuTel 


dAu-I 

2.719 

2.711 

2.680 


dAu-Xe 

2.671 

2.665 

2.642 


dAu-Xe(I) 

2.702 

2.694 

2.654 


dAu-Xe(II) 

2.727 

2.694 

2.684 


dxo ■Xe(I) 

3.337 

3.375 

3.235 


dAu .^i) 

3.711 

3.584 

3.500 


dAu ..^ii) 

3.753 

3.791 

3.633 


dAu..Au{I) 

4.212 

4.163 

4.193 


dAu..Au{II) 

4.143 

4.034 

4.000 


dAu ■Au{III) 

3.998 

3.700 

3.843 


dl-Au-Xe 

160.89 

158.97 

160.26 


dl-Au-Xe(I) 

85.99 

84.27 

86.17 


dl-Au-Xe(II) 

100.78 

100.79 

100.09 


dxe-Au-Xe{I) 

76.77 

78.07 

75.29 


dxe-Au-Xe(II) 

98.02 

98.44 

99.40 


dxe(I)-Au-Xe{II) 

163.69 

170.06 

167.71 


dAu-Xe-Au(I) 

103.23 

101.93 

104.71 


dAu-Xe-Au(II) 

100.25 

97.65 

97.34 


dAu(I)-Xe-Au(II) 

94.86 

86.75 

92.12 

PdTel 


dpd-Xey 

2.636 

2.629 

2.601 


dpd-Xej^ 

2.881 

2.832 

2.844 


dpd-I 

2.762 

2.742 

2.715 


dpd- . . Pd 

4.047 

3.946 

3.952 


dpd- . . Pd 

4.048 

3.910 

3.873 


dxe|| ■ ■ ■ Xey 

3.377 

3.532 

3.385 


dxe ■ ■ ■ Xe 

3.727 

3.772 

3.706 


dp .. I 

3.758 

3.844 

3.809 


dxe|| -Pd-Xey 

79.7 

84.4 

81.2 


dl-Pd-I 

85.74 

89.02 

89.06 


dxc^-Pd-Xej^ 

166.64 

172.64 

168.62 


“Computed by using the lattice parameters given in Tables 1 and 2 of Ref. 15 since in Ref. 15 the bond lengths given in 
Table 3 are not consistent with Tables 1 and 2, as noted by Bahramy et al. (Ref. 5). 

*’Note that in Ref. 15 the bond angles given in Table 3 are consistent with Tables 1 and 2. 
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FIG. SI. (Color online) The plot of the energy difference /S.E = E — {Ems + F'Xe + Ei), cf. Eq. (1), versus the volume for BiTel 
(a), AuTel (b), and PdTel (c). Note that kmax is reduced to 4 for BiTel, and 3 for AuTel and PdTel. It is discernible that a 
good BM fit is achieved with these values of /cmax- 
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FIG. S2. (Color online) The calculated and experimental pressure-volume curves for BiTel. The PBE-calculated and dispersion- 
corrected (PBE-D2) curves are in blue and red, respectively. The empty circles and triangles represent the experimental data 
provided by X. Xi (Ref. 3) and Y. Chen (Ref. 14), respectively. The filled square marks the experimental value of the equilibrium 
volume measured by Shevelkov et al. (Ref. 15). 
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FIG. S3. (Color online) The be planar average of the gradient of the BiTel crystal potential as a function of distance x along 
the a-axis, obtained via semilocal (PBE) and dispersion-corrected (PBE-D2) calcnlations. 



